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Abstract 

We consider correlation functions containing the axial current of one light and 
one heavy quark in the static approximation as well as in full QCD, using the lattice 
regularization. Up to one-loop order of perturbation theory, we study the difference 
between the full and the effective theory in the continuum limit. In the full theory 
we find a term non-analytic in 1/m, revealing the asymptotic character of the \/m- 
expansion. In general, deviations from the m — > oo limit turn out to be small and 
are well described by the first non-trivial terms when m is a factor 2-3 above the 
external scale. We also investigate the mass dependence of discretization errors, and 
find that the behaviour of the correlation functions at finite lattice spacing differs 
significantly from that in the continuum limit when the quark mass is large. 
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1 Introduction 



Heavy quark effective theory (HQET) is an important tool for heavy quark 
phenomenology. As a means of decoupling the light degrees of freedom from 
the mass and spin of heavy quarks (with masses that are much larger than 
Aqcd), it is a powerful method to calculate e.g. the heavy hadron spectrum 
and weak decay form factors. For a review of HQET see [1]. An interesting 
object to study is the axial current of one light and one heavy quark, as the 
B-meson decay constant /b is defined by a matrix element of such a current. 

QCD matrix elements have to be evaluated non-perturbatively, which means 
that the lattice regularization has to be used. As only finite volumes can be 
treated numerically, and the volume must be large enough to accommodate 
the physical system such that finite size effects are negligible, current com- 
puting power enforces a lower bound on accessible lattice spacings, i.e. an 
upper bound on the cutoff. As, on today's computers, this cutoff cannot be 
chosen much larger than the b quark mass, the simulation of bottom fiavoured 
particles is not possible in full QCD, and the use of effective theories like 
HQET becomes a necessity. 

In a previous paper [2], we presented the outline of a method for the non- 
perturbative renormalization of the heavy-light axial current in the static ap- 
proximation, which corresponds to the leading order in the heavy quark ex- 
pansion, where we used Schrodinger functional techniques. In this framework, 
one considers QCD in a finite volume of linear size L and suitable boundary 
conditions. The length scale L takes over the role played by external mo- 
menta Pfj_ in transition amplitudes in infinite volume. HQET emerges in the 
large mass limit at fixed L. When no other scale is present, as it is the case 
in perturbation theory, one has to consider the z ^ oo limit with z — mL 
and some definition of the heavy quark mass, m. The predictive power of 
HQET comes from the fact that the matching of QCD and HQET does not 
depend on the external states (momenta). It is therefore also independent of 
the boundary conditions imposed (and L). For this reason we may use the 
Schrodinger functional to study the relation between the static approximation 
and full QCD. The advantage is that on the one hand gauge invariant, infrared 
finite correlation functions may be computed in perturbation theory (rather 
easily) and on the other hand its lattice-regulated version is very well suited 
for non-perturbative Monte Carlo simulations [3]. 

In this paper, we investigate the difference between the static approxima- 
tion and full QCD in the continuum limit, for the specific correlation func- 
tions defined in [2], both at tree level and at one-loop order of perturbation 
theory. For the first time (to our knowledge), it is explicitly demonstrated that 
gauge invariant, infrared finite correlation functions are described by HQET 
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in the limit m — > oo. Indeed, we observe that the 1/m-expansion is only an 
asymptotic expansion; small non-analytic terms exp(— ^) are present in the 
full theory. 

A second point we look at in this paper is the mass dependence of discretiza- 
tion errors, after applying Symanzik's improvement programme [4] to actions, 
boundary fields, and the axial current. 

We use the notation of reference [2]. This paper is organized as follows. In 
section 2, the basic definitions of the static quark action, both in the continuum 
and on the lattice, are given. Schrodinger functional boundary conditions and 
the correlation functions we use are also explained in that section. In section 
3, the renormalization schemes we use in the static approximation and in full 
QCD are introduced. The investigation of finite mass effects, i.e. the difference 
between the full and effective theories, is carried out in section 4. Section 5 
is devoted to the study of discretization errors and their dependence on the 
heavy quark mass. Our results are summarized in section 6. An appendix 
contains the relevant one-loop Feynman diagrams. 



2 Correlation functions 

2.1 Static quarks 

Static quarks can be treated as two-component objects. To simplify the no- 
tation, they will be described as four-component spinors here, satisfying the 
constraint 



1(1 + 7o). 



(2.1) 



The static quark action is 





Following [5], we discretize the static quark action as 



(2.3) 



X 
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2.2 The Schro dinger functional 



The Schrodinger functional consists of a finite space-time volume TxLxLxL, 
with specifically chosen boundary conditions. Explicitly, the gauge field is 
chosen to be periodic in space and to satisfy Dirichlet boundary conditions in 
time, 

Ak{x)\,,=o = C(x), Ak{x)\,,=T = [PC']{^) , (2.4) 

where C(x) and C"(x) are fixed boundary fields, and P projects onto the 
gauge invariant content of C"(x) (see [6]). Here we choose C — C — 0, leading 
to the boundary conditions 

U{x,k)\,,=o^l, U{x,k)\,,=T^l (2.5) 
for the lattice gauge field. 

The light quark field is chosen to be periodic up to a phase in the three space 
directions, 

V^(x + Lk) = e^V(^), '^{x + Lk) - i:{x)e-'\ (2.6) 
where 6 is kept as a free parameter. 

As for the gauge field, Dirichlet boundary conditions are imposed on the light 
quark field, 

-P+^i(a:^)Uo=o = a(x), P-ipi{x)\^^=T = P\{^), (2.7) 

and 

z^,(x)P_|^o=o = a(x), V^i(x)P+|^o=r = a'(x), (2.8) 
as well as on the static quark field, 

V'h(2;)Uo=o = Ph(x), V'h(2^)Uo=T = Ph'(x). (2.9) 

Note that no projectors are necessary in the heavy quark case because of 
eq. (2.1). For the same reason, we have P^ipi^lx) = 0. Spatial boundary con- 
ditions do not need to be discussed for the static quarks since they do not 
propagate in space. The quark field boundary conditions given above can be 
applied both in the continuum and on the lattice. 
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LxLxL 



LxLxL 



Fig. 1. Schematic drawing of the correlation functions /a (left) and /i (right). The 
dot in the middle of the left diagram symbolizes the axial current at xq, the other 
dots represent the boundary quark fields. 

The Schrodinger functional action, including boundary and 0(a) improvement 
terms, can be found in [6,7,2]. In this paper, it is understood that all quantities 
are calculated using the 0(a) improved action. 

2.3 Correlation functions 

In the following, we study two different correlation functions. The function /a 
contains the time component of a static-light axial current 



where the indices 1 and 2 label two flavours of relativistic quarks. It is cor- 
related with a suitable combination of "boundary quark fields" to yield the 
diagram illustrated in the left of Figure 1 upon Wick-contraction. In addition 
a boundary-to-boundary correlation /i, shown on the right of the figure is 
needed. 

With the boundary quark fields C (see [2] for their proper definition), the 
correlation functions 



stat 
■0 



(2.10) 



or a light-light axial current 



Ao{x) = V'i(a;)7o75V'2(a;), 



(2.11) 



dVd^z«^*(x)Ch(y)75Ci(z)) 



(2.12) 
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and 



/r = d=^ud=^vdVd^z(C/(u)75Ch'(v)Ch(y)75Ci(z)) (2.13) 

can be introduced as in [2]. In addition, the corresponding functions in full 
QCD, 

Uixo) = -\j dVd^z (^(^)C2(y)75Ci(z)) , (2.14) 



and 



h = / d^udVdVd^^z (Ci'(u)75C2'(v)C2(y)75Ci(z)) (2.15) 

are used. They are discretized as in references [8,9], including 0(a) improve- 
ment terms 



^^stat^^stat ^ acf Vii(^.- + V*)7,-75V'h (2.16) 

and 

acA^o = acA\{dl + ^o) {^/'i75^/'2} , (2.17) 
with the improvement coefficients 



cT = glcT^'^ + 0{gtl cT^'^ = -i- X 1.00(1) (2.18) 



and 



CA = glc^^ + O(^o'), c^K = -0.00567(1)Cf, (2.19) 

where the coefficient c^^*^^'* has been calculated in [10] using NRQCD (a cal- 
culation using a different method can be found in [11]), and the coefficient c^^ 
has been taken from [12]. 

With these correlation functions, the ratios 



X,{go,L/a)= ^^^ XL ^ (2.20) 



and 



Yi{go,z,L/a) = ^-^0^ (2.21) 
v/i 

are defined, where the index I means that the 0(a) improved axial currents 
are used, and z parametrizes the mass dependence as defined in subsection 3.1 
below. In the ratios Xi and li, renormalization constants for the light and 
static quark boundary fields cancel as well as a 1/a divergence due to a heavy 
quark mass counterterm in the static approximation (see [2] for details). 

For our perturbative analysis, we expand the ratios to one-loop order, 

Xi(5o, L/a) = x['\L/a) + x['\L/a)gl + O(^o'), (2-22) 

yi(po, z, L/a) = yI'\z, L/a) + Yi'\z, L/a)gl + 0(^o')- (2-23) 
3 Renormalization and matching 



When taking the continuum limit, renormalization of quark masses, couplings, 
and composite operators is required. Bare lattice quark masses tuq need an 
additive renormalization 



mq = mo — rric (3.1) 

due to the chiral symmetry breaking term in the Wilson fermion action, de- 
fining a reduced quark mass mq. The critical mass rric is determined by the 
requirement that in the case of two mass-degenerate quarks the PCAC mass 
is zero at mq = 0. In perturbation theory, the critical mass can be expanded 
as 



me = 5oW^ + 0(^o')- (3-2) 

Depending on the precise definition of the PCAC mass, rric (and therefore also 

m^^)) depend on the lattice spacing. However, in perturbation theory one may 
take the limit where external momentum scales are small compared to the 
cutoff l/a, and extrapolate the coefficient m^^^ to its continuum limit. This 
gives [12] 
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am(^) = -0.2025565(1) x Cp, 



(3.3) 



a value which will be used also in Sect. 5 at finite values of a. Throughout this 
paper, we use rriq — for the light quark field ipi. 

3.1 Renormalization schemes 

Unlike the case of full QCD, which is known to be renormalizable at all or- 
ders of lattice perturbation theory [13], renormalizability of the static theory 
has not yet been proven. However, a large number of calculations has been 
performed in the static approximation, and none of these has given a reason 
to doubt renormalizability. Throughout this paper, we assume that the static 
theory is renormalizable "as usual" by adding local counterterms to action 
and composite fields. 

At one-loop level, the ratio X is divergent in the continuum limit, i.e. a scale 
dependent renormalization is required. The one-loop coefficient X^^^ is ex- 
pected to diverge logarithmically, and a renormalization scheme, the lattice 
MS scheme can be defined by subtracting this logarithm, 

= ZffMXi{a/L), (3.4) 

ZtfM = I - ^oh^MgL, (3.5) 
where 

gL-g'o + oig'o), (3.6) 

and fi is the renormalization scale. The one-loop anomalous dimension of the 
static axial current 



To = -i^ (3.7) 

does not depend on the renormalization scheme and is known from [14,15]. 

The ratio Yi{z,a/L) remains finite in the continuum limit as long as the 
coupling and the heavy quark mass are made finite by renormalization. The 
renormalization condition used here is that, in the continuum limit, they equal 
the MS coupling and quark masses respectively, at the renormalization scale 
1/L. At one-loop level, this means that the renormalized quantities 
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9ws = 9l + 0{gt) 

and 



(3.8) 



^2,MS = ^m,Ms(l + a^mmq,2)mq,2 (3.9) 

are used, where mq^2 is the reduced quark mass of the second quark flavour 
(remember that the mass of the flrst flavour is chosen to be zero) , 



Z^,m = 1 - [^H^/L) - 0.122282CF ) gi^s + O(^Ms) (3-10) 
is known from [16], and the improvement coefficient 



fcm = - 0.07217(2)Cf^2 + O(^o') (3-11) 
has been calculated in [9]. 

Taking the continuum limit means to send the lattice spacing to zero while 
keeping the physical box size fixed. One must thus take the limit £ — > 0, 
while g^ig and z = Lm2^ are kept at fixed values. 

A finite renormalization of the light-light axial current is possible, and we 
define a renormalized current (Aca)o by requiring that it satisfies the current 
algebra relations in the case ;z = (see [17,8]). This renormalization condition, 
together with the MS definition of the coupling and the quark mass, gives a 
renormalized ratio 



Yc^{z, a/L) = Za(1 + \ahp^m^,2)Yi, (3.12) 
with the perturbative renormalization constant [16] 

Za = 1 - zi^C^gl + O(^o'), Z^K = -0.0873435, (3.13) 
and the improvement coefficient [9] 

6a = 1 + 0.11414(4)Cf^o' + 0(^o')- (3-14) 

In principle, IcA should be considered also as a function of g"^. Since at l-loop 
order the dependence on the coupling is rather trivial, we have not indicated 
it explicitly. 
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3.2 Matching of full and effective theories 



By construction, the effective theory is expected to describe the physics of full 
QCD in the limit m — > oo, when this limit is taken appropriately [18]. Beyond 
the classical level this physics boundary condition translates into a matching 
condition for the two theories, fixing a renormalization scheme ("match") for 
the effective theory. The details of the matching conditions should of course 
be irrelevant, as long as they fix the renormalization constants. To match the 
axial current, we choose the condition 

^match(m2,Ms) = Yca{z, a/L) + 0(1/ z) + 0((a/L)2), (3.15) 
which means that the effective theory and full QCD are matched at the scale 

The renormalized static axial current in this scheme is related to the current 
in the lattice MS scheme by finite renormalization, i. e. 

^match(A^) = (xr)^iat(A^) + 0{{a/Lf), (3.16) 

with 

ixr) = l + BrgL + 0{gL), (3.17) 
where B^^^^ is a finite constant. Using 

= iUL) - 70 ln(^)X(°) (a/L) + 0((a/L)^) , (3.18) 
this constant can be calculated from the matching condition eq. (3.15), 

Bt^X('^ {a/L) = F« {z, a/L) - X« ( 1/L) + 70 Hz)X^'^ {a/L) 

+0{l/z) + 0{{a/L)^). (3.19) 

Here it is understood that all quantities are extrapolated to the continuum 
limit first, and then the extrapolation to I/2; = is performed to determine 
^stat practical terms, we define a quantity B^^^{z) as the continuum limit 
of 



Bt\z, a/L) =70 In(^) + ^mj^) l^i^'^^^' "/^^ 
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Fig. 2. Bf^^{z,a/L) at 6* = 0.5, for 2 = 8 (squares) and z = 12 (triangles). The 
bigger symbols are the data for Lla = 32. 



+4'VW(z,a/L)-X«(a/L)}. 



(3.20) 



The continuum limit of X*^^^ and Y^^"^ is calculated separately by fitting the 
numerical data with a function of the expected form [19] 



+ /i3^ + /i4^1n(l/0 + O(l//^), 



/ — L/a, 



(3.21) 



using the fit procedure described in reference [20]. This reference also explains 
a method to estimate the errors of the fit parameters, which we used in our 
calculations. For d — 0.5, we used eq. (3.21) without the 0{{a/Vf) terms for 
z > b and included those terms in all other fits. For some parameters, the 13 
values resulting from our numerical data and the fits are shown in Fig. 2. 

It is expected that 



stat 
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Fig. 3. Calculation of the matching coefficient B^^* by extrapolation to 1/z = 
for 9 = 0.0 (triangles) and 6 = 0.5 (squares). Where not indicated the errors are 
smaller than the symbol sizes. 



z 



1 

z 



o(- 



(3.22) 



where 0(1/2;^) means terms of the form 1/z'^ and ln{z)/z'^, and the constants 
/i and /2 can be determined from a fit to the numerical data. Figure 3 shows 
the data for a set of values of z between 1 and 12 at 6' = 0.0 and between 
1 and 16 at ^ = 0.5, as well as fits of the form (3.22), where only the filled 
symbols have been included in the fits. The results for the fit constants are 



e 


Dstat 


/i 


h 


0.0 
0.5 


-0.136(3) 
-0.137(1) 


-0.05(2) 
0.011(5) 


-0.04(2) 
-0.054(5) 



which means that our final result for Bf^^ is 
^stat ^ _o.i37(l). 



(3.23) 
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This constant is already known from [21], where it has been calculated from 
matrix elements between quark states, using a gluon mass as infrared regula- 
tor. Extracting the relevant terms from that calculation, one gets 



(3.24) 



in agreement with our result. 

It is more illuminating to turn this logic around: If we take the matching con- 
stant S^^* from matrix elements between quark states, the correlation func- 
tions calculated from the static effective theory are the infinite mass limits of 
the corresponding correlation functions in full QCD, and the mass-dependence 
of the QCD correlation functions at one loop level has the form expected from 
HQET. To our knowledge this is the first explicit demonstration of this fact 
for gauge invariant, infrared finite correlation functions. 



4 Heavy quark expansion as an asymptotic series — finite mass 



As Xmatch IS dcfiucd by matching to IcA, the static theory and QCD with 
relativistic quarks can be compared directly, when these two renormalization 
schemes are used. To calculate the actual size of the finite mass effects, i.e. 
the difference between the ratio Yqa and the static limit, we define 



At tree level, the continuum limits of X^°^ and Y^^^ can be calculated ana- 
lytically. 



effects 



A(^) 



^Ca(^) — -^match (^2,Ms) 
^match (^^2,1^13) 



(4.1) 



and expand 



A(z)=A(o)(z) + AW(%|g + 0(^y. 



(4.2) 




2Ri 



3uJi 



{l + e-n + Oiia/Lf), 



(4.3) 



y^°)(z,a/L) 



1 



{{z + u;2) + iz- u;2)e-^'^}X^^\a/L) 
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Fig. 4. Y{z) at tree level (solid curve), for T = L and two values of 9, with the 
expansion up to order 1/z (long dashes), order (l/z)'^ (short dashes), and order 
{l/zf (dotted curve . 



+ , , „ „ (1 - e-'^^^)(l - e-'^^^) + 0((a/L)2), (4.4) 



where 



Ri^ui{l + e-'^^^), R2 ^ z + u;2 - {z - u;2)e-^^'\ (4.6) 

Eq. (4.4) shows that the continuum hmit of Y^^^ is non-analytic in 1/z, which 
means that an expansion in 1/z in the sense of a convergent Taylor series is 
not possible, the expansion is at most asymptotic. However, the non-analytic 
terms in eq. (4.4) are suppressed by terms of the form e~'^'^'^, and in the static 
limit z — > 00, the ratio Y^^^ equals the ratio X^^^ calculated from the static 
quark action eq. (2.2). While y(°)(z) = X^^^ for ^ = 0, there are finite mass 
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effects for other values of 0. For 9 — 0.5 and 6 = 1.0 tlie resulting curves for 
Y^^^ as a function of I/2; are shown in Fig. 4, together with the expansion up 
to order {1/zy. The I/2; — > limit of these plots gives the respective value 
of X^^\ It is clearly visible that the heavy quark expansion breaks down for 
z < 2, which corresponds to m2^MS < 2/iv. 

From eq. (4.3) and eq. (4.4), we get the tree level quantity A^^\ The result at 
T = L is that the finite mass effects for 2; > 1 are in the range of 0-14 per 
cent for 6 = 0.5 and 0-40 per cent for 6 = 1.0. Details can be read off from 
Fig. 4. 

At one-loop level, the correlation functions are evaluated numerically as de- 
scribed above, and A^^^z) is calculated from the numerical data shown in 
Fig. 3. For the set of parameters used in Fig. 3, A*^^) lies between and 5 per 
cent for 9 = 0.0, and between and 2 per cent for 9 = 0.5. 
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Fig. 5. Discretization errors in the ratio Y{z) at tree level. The circles are z = 3 
data, the triangles z = A, the squares z = 8, and the pentagons z = 16. 
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Fig. 6. Discretization errors in the ratio Y{z) at one-loop level. The circles are z = 3 
data, the triangles z = A, the squares z = 8, and the pentagons z = 16. 



5 Cutoff effects 



As we use an 0(a) improved action, improved boundary fields, and improved 
operators in our lattice calculation, we expect discretization errors of O(a^). 
These are expected to depend on the heavy quark mass. To shed some light 
on the mass dependence of the discretization errors, we define 



5{z, am^ j^) = 
and expand again. 



YcA{z,a/L)-YcA{z,0) 
Yca{z,Q) 



(5.1) 



5 = 5(°)+c5«^,|-3 + 0(^y. (5.2) 

We calculate 5*^°^ and 5^-^^ for a range of different values of z and a/L, and print 
them as a function of (a^rigMs)^ in Fig. 5 and Fig. 6. The figures show that (for 
our choice of parameters) the discretization errors are roughly proportional 
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Fig. 7. The matching coefficient B^^* in the continuum hmit (squares) and at 
L/a = 32 (triangles) for 6 = 0.5. The dotted line is the fit explained in Sect. 3. 



to (0^2 ^13)^. This is a significant problem in heavy quark extrapolations at 
finite lattice spacing, as due to the mass dependence of the discretization 
errors, there is a certain risk that the slope being used in the extrapolation 
contains a sizeable contribution from discretization errors. To illustrate this a 
bit further. Fig. 7 shows the quantity 5^^* introduced in Sect. 3, both in the 
continuum with the fit described above, and at finite lattice spacing L/a = 32. 
Here again it can be seen that the discretization errors increase with the value 
of the heavy quark mass, and for z > 8 (which means am2y[g > 1/4)) the 
behaviour of the curve at finite lattice spacing changes dramatically. That 
this is really due to a cutoff effect can be seen from Fig. 2, where the data 
points for L/a = 32 have been marked by slightly bigger symbols. This shows 
that 0(a) improvement breaks down for too large values of arn2 ms' heavy 
quark extrapolations should be performed in the continuum limit rather than 
at finite lattice spacing. A detailed explanation of the improvement breakdown 
for heavy quarks can be found in section 4 of ref. [22], where this effect is 
discussed in the context of the Schrodinger functional renormalization of the 
QCD gauge coupling. 
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A comparison between the values of 5(^,ar?T2Ms) for Ms < 1/2 and A(z) 
for 2; > 2, both at tree and one-foop level, shows that the discretization errors 
and the finite mass effects are of comparable size, i.e. of the order of a few per 
cent. 



6 Conclusions 



The main point of this paper is a comparison between the axial current in 
QCD with finite quark masses, and in the static approximation at one-loop 
order of perturbation theory. To achieve this, we have defined gauge invariant, 
infrared finite, finite- volume correlation functions with Schrodinger functional 
boundary conditions, and constructed ratios of these correlation functions to 
cancel renormalization constants except for that of the axial current. Further- 
more, the divergence due to the residual mass term in the static theory cancels 
in our ratios. 

Using these ratios, and assuming renormalizability of the static theory, we 
have defined a renormalization scheme for the static-light axial current by 
matching to the full theory, such that the renormalized static ratio equals the 
heavy-quark limit of the corresponding ratio calculated from full QCD. This 
renormalization scheme allows a comparison between the full and the static 
theory. Taking the matching constant 5^^* from one specific matching condi- 
tion, we have verified explicitly that the heavy quark limit of an independent 
correlation function is described by HQET. For quark masses down to l/L, 
where L is the extent of the space-time box, the finite-mass effects in our ratio 
are in the range of a few per cent, both at tree level and at one-loop order of 
perturbation theory. 

As we performed our calculations in the lattice regularization, we were also 
able to investigate discretization errors in our ratio of correlation functions, 
both at tree and one-loop level. We have shown that (for our choice of pa- 
rameters), the tree level and one-loop coefficients approach the continuum 
limit at a rate roughly proportional to (am)^, if 0(a) improvement is applied 
to the action, the boundary fields, and to composite operators. We observe a 
significant difference between continuum limit data and results at finite lattice 
spacing, with a qualitative difference in the behaviour of the two curves for 
am > 1/4. This demonstrates that extrapolations in the heavy quark mass 
should be carried out in the continuum limit rather than at finite lattice 
spacings. 
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A Feynman diagrams 



The expansion of the correlation functions /a and /i is described in detail 
in ref. [12]. Following that method, we get the one-loop Feynman diagrams 
shown in Fig. A.l and Fig. A. 2. 
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1a 
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1b 
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Fig. A.l. One loop diagrams contributing to /a (T/2). 



19 



T-a T 
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Fig. A. 2. One loop diagrams contributing to fi. The dotted lines denote the link 
from T — a to T. 



References 

[1] M. Neubert, Phys. Kept. 245 (1994) 259, hep-ph/9306320, and references 
therein. 

[2] ALPHA, M. Kurth and R. Sommer, Nucl. Phys. B597 (2001) 488, hep- 
lat/0007002. 

[3] K. Jansen et al., Phys. Lett. B372 (1996) 275, hep-lat/9512009. 

[4] K. Symanzik, Nucl. Phys. B226 (1983) 187. 

[5] E. Eichten and B. Hill, Phys. Lett. B234 (1990) 511. 

[6] M. Liischer, R. Narayanan, P. Weisz and U. Wolff, Nucl. Phys. B384 (1992) 
168, hep-lat/9207009. 

[7] S. Sint, Nucl. Phys. B421 (1994) 135, hep-lat/9312079. 

[8] M. Liischer, S. Sint, R. Sommer and H. Wittig, Nucl. Phys. B491 (1997) 344, 
hep-lat/9611015. 

[9] S. Sint and P. Weisz, Nucl. Phys. B502 (1997) 251, hep-lat/9704001. 



20 



[10] C. Morningstar and J. Shigemitsu, Phys. Rev. D57 (1998), hep-lat/9712015. 

[11] K.I. Ishikawa, T. Onogi and N. Yamada, Nucl. Phys. Proc. Suppl. 83 (2000) 
301, hep-lat/9909159. 

[12] M. Liischer and P. Weisz, Nucl. Phys. B479 (1996) 429, hep-lat/9606016. 

[13] T. Reisz, Nucl. Phys. B318 (1989) 417. 

[14] M.A. Shifman and M.B. Voloshin, Sov. J. Nucl. Phys. 45 (1987) 292. 

[15] H.D. Politzer and M.B. Wise, Phys. Lett. B206 (1988) 681. 

[16] E. Gabrielh, G. Martinelh, C. Pittori, G. Heatlie and C.T. Sachrajda, Nucl. 
Phys. B362 (1991) 475. 

[17] M. Bochicchio, L. Maiani, G. Martinelh, G.C. Rossi and M. Testa, Nucl. Phys. 
B262 (1985) 331. 

[18] T. Mannel, W. Roberts and Z. Ryzak, Nucl. Phys. B368 (1992) 204. 
[19] M. Liischer and P. Weisz, Nucl. Phys. B266 (1986) 309. 

[20] ALPHA, A. Bode, P. Weisz and U. Wolff, Nucl. Phys. B576 (2000) 517, hep- 
lat/9911018. 

[21] A. Borrelh and C. Pittori, Nucl. Phys. B385 (1992) 502. 

[22] S. Sint and R. Sommer, Nucl. Phys. B465 (1996) 71, hep-lat/9508012. 



21 



